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Abstract. Three-dimensional magnetohydrodynamical large eddy simula- 
tions of solar surface convection using realistic model physics is conducted. The 
effects of magnetic fields on thermal structure of convective motions into radia- 
tive layers, the range of convection cell sizes and penetration depths of convection 
is investigated. We simulate a some portion of the solar photosphere and the 
upper layers of the convection zone, a region extending 30 x 30 Mm horizontally 
from Mm down to 18 Mm below the visible surface. We solve equations of 
the fully compressible radiation magnetohydrodynamics with dynamical viscos- 
ity and gravity. For numerical simulation we use: l)realistic initial model of Sun 
and equation of state and opacities of stellar matter, 2) high order conservative 
TVD scheme for solution magnetohydrodynamics, 3) diffusion approximation for 
solution radiative transfer 4) calculation dynamical viscosity from subgrid scale 
modelling. Simulations are conducted on horizontal uniform grid of 320 x 320 
and with 144 nonuniformly spaced vertical grid points on the 128 processors of 
supercomputer MBC-1500 with distributed memory multiprocessors in Russian 
Academy of Sciences. 



1. Introduction 

Convection near solar surface has strongly non-local and dynamical charac- 
ter. Hence numerical simulation provide useful information on the structure 
spatial scales by convection and help to construct consistent models of the 
physical processes underlying the observed solar phenomena. We investigate 
effects of compressibility and weak of magnetic field on formation non-local 
structure of convection using realistic physics and conservative TVD numeri- 
cal scheme of Godunov type. The previous simulations were confined by small 
computational domain an d studied processes on scales order size of granulation 
[Stein Nordlund I ((2006)] . In order to investigate collective interaction of con- 
vective modes different scales and process of formation of supergranulation we 
conducted calculation in three dimensional computational box by size 30 Mm in 
horizontal direction and by size 18 Mm in vertical direction. 



2. Numerical method 

We take distributio n of the main thermodynamic variables by radius due to Stan- 
dard Solar Model [Christensen-Dalsgaardl (J2003J)] with parameters (X, Z, a) = 
(0.7385,0.0181,2.02), where X and Y are hydrogen and helium abundance by 
mass, and a is the ratio of mixing length to pressure scale height in convection 
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region . We use OPAL opacities and equations of state for solar matter [iglesiasl 
(199(1)]. 



We solve fully compressible nonideal magnetohydrodynamics equations: 
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where E = e + pv 2 /2 + B 2 /8tt is the total energy, Q ra d is the energy trans- 
ferred by radiation and r is viscous stress tensor. 

We assume that small scales are independent of resolved scales (Large Eddy 
Simul ations) and rate dissi pation is defined from buoyancy and shear production 
terms |Canuto et alJ (|l994l )]. The numerical method that we us e d was an explicit 
Godunov-type conservative TVD difference scheme |Yee et all (1990)] 
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Flux along each direction, for example x, was defined by local-characteristic 
method as follows 

Fi+l/2,j,k = 9 [Fij,k + Fi+i,j,k + ^1+1/2^+1/2] 



where -Ri+1/2 is matrix whose columns are right eigenvectors of dF/dU 
evaluated at generalized Roe average for real gases of Uij^ and C/j+ij^. The 
Wj+i/2 is the matrix of numerical dissipation. Term Sij t k is accounted effect of 
gravitation forces and radiation. 

The one step of time integration is defined by Runge-Kutta method Shu &i Qsherl 
(1988)] as 
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The scheme is second order by space and time. For approximation vis- 
cous terms we used central differences. For evaluation radiative term in energy 
equation we used the diffusion approximation 
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We use uniform grid in x-y directions and nonuniform grid in vertical z one. 
We apply periodic boundary condition in horizontal planes and choose on the 
top and bottom as follows 

dp/dz = pg z ,P = p(p), e = const 



By = 0,dB z /dz = 



In initial moment magnetic field equal 50 G and has just one vertical compo- 
nent. For the magnetic diffusitvity we take constant value r\ = 1.1 x 10 11 cm 2 sec _1 



3. Results 

On the figures 1-4 results of development convection after 12 solar hr MHD nu- 
merical simulation are shown. We founded that magnetic field concentrate on 
the boundary of convective cells in forms magnetic flux and sheets. Diverging 
of convective flows from centre supergranular expell weak magnetic field on the 
edges of convective cell. Average size of supergranular celss is 10-20Mm with 
lifetime about 8-10 hours. We aplly procedure averaging by interval time two 
hour and find value of velocity from centre supergranular equal to 1-1.5 km/sec. 
In places action of strong magnetic field with strength 700-900 G we observe 
effect of suppression of convection and decreasing of fluctuation of temperature. 
Magnetic pressure in regions concentration of magnetic flux prevent inflow of 
matter. Transfer of radiation energy in these places is suppressed. We founded 
from simulation that maximum of value of magnetic field in computational do- 
main equal to 1300 G. 

Inside of supergranular we have usual picture of evolution of convection on 
scale of granulation with average sizes of cells about l-2Mm and lifetimes about 
10 minutes. Here we see wider upflows of warm, low density, and entropy neu- 
tral matter and downflows of cold, converging into filamentary structures, dense 
material. We observe continuous picture formation and destruction of gran- 
ules. Granules with highest pressure grow and push matter against neighboring 
granules, that then shrink and disappear. Ascending flow increases pressure in 
center of granule and upflowing fluids decelerates motion. This process reduce 
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heat transport to surface and allow material above the granule to cool, become 
denser, and by action gravity to move down. We observe formation new cold 
inter granule lane splitting the original granule. 

From figure 4 we see distinctly existence three different regions development 
of convection. In near of solar surface to the depth 4 Mm we founded zone 
of turbulent convection. In this region cold blobs of matter move down with 
velocity in maximum about 4 km/sec with maximum Mach number equal 1. 
The downdrafts has different and very complicate vertical structure. Some ones 
travel small distance from surface and become weak enough to be broken up by 
the surrounding fluid motion. Other ones conserve motion with high velocity and 
move on distance about 6 Mm. We observe that different nature such behavior 
is due to initial condition of formation downdrafts. In place of confluence of 
convective cells to one point more energy is released. 

In region from 5Mm to 8Mm of depth we reveal more quiet character of 
convective flow than in turbulent zone. Below 8Mm we see clear separate large 
scale density fluctuations and streaming flow of matter similar jets with average 
velocity about lkm/sec. The magnetic field in these places has values about 
300G. Distance beetween different narrow jets gives size of supergranular cells. 
Places of intersections of path jets with horizontal plane are vertices of huge 
convective cells. From distribution iso-surface of magnetic field with strength 
500 G we have discovered that pumping of magnetic field occupy more part of 
all computational domain up to bottom boundary. 

Magnetic field near solar surface in our numerical simulation has very com- 
plicate structure. Inside supergranular cells we see formation, growth and evo- 
lution with time many arising to surface loops of magnetic field. There is places 
in vertices of huge cells with vorticity motion that provide quick rate magnetic 
helicity transport across solar photosphere(Fig.5,6). On the boundary super- 
granular cell magnetic field has just vertical component. In these parts we 
observe quick changing sign and big values of current helicity (Fig. 7). 
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Figure 2. Image of contours of the vertical component of magnetic field at 
the same plane as that for Figure 1. 



Ustyugov 




Figure 4. Image of contours of the fluctuations density and field of velocity 
in vertical plane. 
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Figure 6. Image of contours of the magnetic hclicity transport rate and field 
of velocity in small part computational domain at the same plane as that for 
Figure 5. 
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Figure 7. Image of contours of the current helicity at the same plane as 
that for Figure 1. 



